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We use the POWHEG formalism in the Herwig++ event generator to match QCD 
real-emission matrix elements with the parton shower for a range of decays relevant 
to Beyond the Standard Model physics searches. Applying this correction affects the 
shapes of experimental observables, changing the number of events passing selection 
cuts. To illustrate this effect we consider the invariant-mass distribution of dijets 
produced in the decay of the lightest Randall-Sundrum graviton and jet measure and 
transverse momentum distributions for decays in Supersymmetry. We consider only 
the effect the POWHEG correction has on the shapes of distributions and ignore the 
global normalisation factor required to correct the total widths and branching ratios 
to next-to-leading order accuracy. 
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1 Introduction 

For Beyond the Standard Model (BSM) scenarios with additional new particles, 
the decays of these particles determine the experimental signals we would observe 
at the LHC If the new particles have a well separated mass spectrum, long 
decay chains will occur when a heavy new particle is produced. For decays 
involving coloured particles, hard QCD radiation at each step in the decay chain 
will alter the structure of the event and therefore the number of events passing 
experimental cuts. The effects of radiation are also important in models with 
degenerate new particle mass spectra where decay chains are typically limited to 
one step. Searches for these compressed spectra scenarios look for events where 
hard radiation in the initial-state shower recoils against missing transverse energy 
in the final state to give an observable signal [I]. The emission of hard QCD 
radiation in the decay of new particles could either enhance or reduce this effect 
and so must be taken into account. Hence, accurate simulation of hard radiation 
in the decays of BSM particles is necessary in order to optimise searches for new 
physics. 

Monte Carlo event generators use fixed-order matrix elements combined with 
parton showers and hardonization models to simulate events. In the Herwig++ 
event generator |2]|3) the decays of unstable fundamental particles are treated sep- 
arately from the hard process which produced them, prior to the parton shower 
phase, using the narrow width approximation. Decays are generated using the 
algorithm described in El which ensures spin correlations are correctly treated. 
The parton shower utilises an approximation that resums the leading collinear and 
leading-colour soft logarithms |5| and so does not accurately describe QCD radi- 
ation in the regions of phase space where the transverse momenta of the emitted 
partons are high. The Positive Weight Hardest Emission Generator (POWHEG) 
formalism [6] is one method that allows the simulation of high transverse momen- 
tum (hard) radiation to be improved by using the real-emission matrix element 
to produce the hardest emission in the shower. This approach affects both the 
overall cross sections for inclusive processes and results in non-global changes 
to the shapes of distributions sensitive to the hardest emission. In particular, 
non-global changes to observables such as jet transverse momenta are important 
since they can impact on the proportion of events passing selection cuts in new 
physics searches. Since BSM signals often consist of only a few events this can 
in turn result in significant changes to the exclusion bounds that can be set. 

The POWHEG formalism has been successfully applied to a wide range of 
hard production processes, for example [7f|27j, and particle decays 28,29 in 



the Standard Model (SM) as well as selected BSM processes 30-34 . Next-to- 



leading order (NLO) corrections to BSM particle decays have also previously been 



studied, for example in 35 where the super-QCD correction to the decay q qx 
was calculated. In this work we present results from the implementation of the 
POWHEG method in Herwig++ for a range of decays relevant for new physics 
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searches. A similar approach based on generic spin structures is used to apply a 



matrix-element correction to hard radiation in particle decays in PYTHIA 36 . 

The POWHEG formalism will be reviewed in Sect. [2] and in Sect. [3] details of 
the decay modes implemented will be given. Results for the decay of the light- 



est graviton in the Randall-Sundrum (RS) model p7] will be given in Sect. 4.1 



Finally, results from a selection of decays in the Constrained Minimal Supersym- 



metric Standard Model (CMSSM) will be presented in Sect. 4.2 



2 POWHEG Method 



In this section a brief outline of the POWHEG method is given. Further details 
can be found in 



38 



In the conventional parton shower approach the inclusive differential cross 
section for the highest transverse momentum emission from an N body process 
is given by, 

da PS = B(^ N )d^ N [A(p Tmin ) + A{p T )d$ R V) . (1) 

Here we are considering a parton shower ordered in terms of the transverse mo- 
mentum of the additional parton, p?- B is the Born-level matrix element squared 
including the relevant flux factor, $tv are the iV-body phase-space variables of 
the leading-order (LO) process and Q R are the phase-space variables of the addi- 
tional parton. V is the unregularized Altarelli-Parisi splitting kernel and PTmin is 
a transverse momentum cut-off introduced to regularise the infra-red (IR) diver- 
gences in the splitting function. The probability of evolving between the scales 
PTmax and pt without emitting a parton is given by the Sudakov form factor 



A (pa 



exp 



Pro 



d$ R V (pr - PTn 



PT 



(2) 



The square bracket in Eq. [T] integrates to unity which ensures that the total cross 
section is given by the LO result. 

In the POWHEG approach the inclusive differential cross section for the hard- 
est emission is given by the QCD NLO differential cross section, that is, 



da PO = B($ N )d$ 



N 



A PO (p Tm - m ) + A™(p T )d$ R 



PO/ 



N 



(3) 



where B(§ N ) is defined by, 



V($ N )+ / C($ N ,$ R )d<5> R 



+ / [R($ N , $ R )d$ R - C($ N , $ R )d<$> R ] 



(4) 
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The real emission contribution, R(<&n,<&r), corresponds to the radiation of an 
additional parton from the LO interaction and the virtual contribution, V(<&n), 
comes from the 1-loop correction to the LO process. C($>n,Qr) is a counter 
term with the same singular behaviour as the real and virtual contributions and 
is introduced to ensure the two square brackets in Eq. [4] are separately finite. 
Typically the counter term is rewritten as a sum of dipole functions, 2^, each 
of which describes the singular behaviour of the matrix elements in a particular 
region of phase space. The Sudakov form factor appearing in Eq. [3] is, 

A PO (p T ) = exp (- jT maX d* R $§0) 6 ( PT - p Tmm )) . (5) 

As with the conventional parton shower approach, the square bracket in Eq. [3] 
will integrate to unity and hence the total inclusive cross section will be given by 
the NLO result. The real-emission matrix element in Eq. [5] can be rewitten using 
dipole functions, 

R( $ N ,$ R )Vi 



fl($*,**) = 2j v w;r . ( 6 ) 



such that Eq. [5] becomes a product of Sudakov form factors, each of which de- 
scribes the non-emission probability in a particular region of phase space specified 
by the dipole function TV 

When applying the POWHEG method to a parton shower ordered in trans- 
verse momentum the hardest emission is generated first using the POWHEG 
Sudakov form factor in Eq. [5j Subsequent emissions are generated with the nor- 
mal parton shower Sudakov given in Eq. [2j with the requirement that no emis- 
sions generated using this form factor have higher transverse momentum than 
the hardest emission coming from However in Herwig++ the shower 

is ordered using an angular evolution variable (2,39 rather than the transverse 



momentum and as a consequence the hardest emission may not be generated 
first. The POWHEG approach can be reconciled with angularly ordered showers 
by splitting the shower into several steps (6l. First the hardest emission in the 
shower is generated using the POWHEG Sudakov form factor and the value of 
the evolution variable corresponding to this emission is determined. A truncated 
parton shower is then used to generate emissions for values of the evolution vari- 
able above this scale. The hardest emission is inserted and the shower is allowed 
to continue down to the IR cut-off of the evolution variable. Finally, any emis- 
sions generated by the shower that are harder than the one generated using the 
POWHEG Sudakov form factor are vetoed. 



3 Decays of BSM Particles 



5 



3 Decays of BSM Particles 

As discussed in Sect. [TJ it is important that the simulation of QCD radiation in 
the decays of BSM particles is done in the most accurate way possible. In this 
work we present results illustrating the effect of consistently matching the QCD 
real-emission matrix element with the parton shower in Herwig++ through the 
POWHEG formalism. This technique has been applied to a range of decays that 
occur in most of the well studied BSM scenarios. Table Q]shows the combinations 
of incoming and outgoing spins for which this method is used and each spin 
structure is implemented for the colour flows given in Table [2] However, models 
with coloured tensor particles are beyond the scope of this work and therefore 
decays involving incoming tensor particles were limited to colour flows in which 
the tensor is a colour singlet. 

The LO and real-emission matrix elements appearing in the POWHEG Su- 
dakov form factor in Eq. [5] are calculated using helicity amplitude methods to 
correctly incorporate spin correlations jl]. The dipole functions, T>i, that appear 



in Eq. |6|are defined as in the Catani-Seymour dipole subtraction method 40,41 
In this approach dipoles describing quasi-collinear radiation from massive vector 
bosons are not well defined. Hence the Fermion-Fermion- Vector, Scalar-Scalar- 
Vector and Tensor- Vector- Vector decays are limited to the situation where any 
final-state coloured vector particles are massless. The Vector-Fermion-Fermion 
and Vector-Scalar-Scalar decays do however include radiation from massive in- 
coming vector particles. Decays are performed in the rest frame of the decaying 
particle [2] and therefore the dipole describing the singular behaviour of this 
particle will only contain a universal soft contribution. This is a well defined, 
spin- independent function. 

Finally, in this work we focus solely on the effect of the POWHEG correction 
on the shapes of observable distributions and have not implemented the global 
normalisation factor coming from the presence of B rather than B in Eq. [3j 
In most cases the partial widths and branching ratios used in Herwig++ are 
calculated by an external program and so the implementation of the B function 
will not affect the results. Even in cases where the calculation of the widths 
and branching ratios is performed in Herwig++ the impact of the POWHEG 
correction on the real radiation will be the dominant effect. 



4 Results 

4.1 Randall-Sundrum Graviton 

The effect of applying the POWHEG correction to the decay of the lightest RS 
graviton was investigated using the Herwig++ implementation of the RS model. 
LHC proton-proton collisions with a centre-of-mass (CM) energy of = 8 TeV 
were simulated. The lightest graviton, G, was produced as a resonance and 
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Tab. 1: Spin combinations for which 
the POWHEG correction 
has been applied. 
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Tab. 2: Colour flows for which the 
POWHEG correction has 
been applied. 



allowed to decay via G —> gg and G —> qq for q = u,d,s,c,b. The mass of 
the graviton was chosen to be mg = 2.23 TeV which corresponds to the lower 
bound on the allowed graviton mass for the coupling k/M p i = 0.1 in 42 . An 
analysis based on the ATLAS experiment's search for new phenomena in dijet 
distributions 43 was then carried out. Jets were constructed using the FastJet 



44 implementation of the anti-fcf algorithm 45 with the energy recombination 

> 4.4 were discarded, 
Events with less than 



scheme and a distance parameter R = 0.6. Jets with \y 
where y is the rapidity of the jet in the pp CM frame, 
two jets passing this constraint were vetoed. The rapidities of the two highest 
Pt jets in the pp CM frame are given by y\ and y?,. In the dijet CM frame 
formed by the two hardest jets their corresponding rapidities are y* and — 
where = \{y\ — 1/2)- Events not satisfying < 0.6 and 1 2/1,2 1 < 2.8 were 
discarded. The dijet invariant mass, rrijj, was formed from the vector sum of the 
two hardest jet momenta and events were vetoed if rrijj < 1.0 TeV. 

The dijet mass distribution after the above selection cuts were applied is shown 
in the left hand plot in Fig. [T] The blue dashed line shows the invariant mass 
distribution for the LO matrix element combined with the parton shower while 
the black line shows the result with the POWHEG correction to the graviton de- 
cay applied. Including the POWHEG correction causes a significant decrease in 
the number of events in the region 2.1 TeV < rrijj < 2.3 TeV. This effect is high- 
lighted in the right hand plot in Fig. [T] which shows the dijet mass distribution 
in this range. The number of events in this region decreases by approximately a 
factor of 1.6 when the POWHEG correction is applied. Of the 10 7 events gener- 
ated 2.4% pass the event selection cuts for the standard parton shower approach 
whereas only 1.5% pass when the POWHEG correction is applied. In the con- 
ventional parton shower approach, the majority of the graviton's momentum will 
be carried by the two partonic decay products. When the POWHEG correction 
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Fig. 1: Dijet invariant mass distribution for the lightest RS graviton decaying 
to jets. The left-hand plot shows the distribution in the full range 
1.0 TeV < rrijj < 2.5 TeV while the right-hand plot emphasises the effect 
on the peak region 2.1 TeV < rrijj < 2.3 TeV. The mass of the graviton 
was vhq = 2.23 TeV and the coupling k/M p i = 0.1. LHC events were 
simulated with yfs = 8 TeV. 



is applied, the highest pt emission in the shower will typically be quite hard 
and so a significant fraction of the the graviton's momentum will be missed by 
considering the invariant mass of only the hardest two jets, therefore shifting the 
distribution to lower values of m,-,-. 



4.2 Constrained Minimal Supersym metric Standard Model 



In addition to the results presented in Sect. |4.1| the effect of the POWHEG cor- 
rection was also studied in the context of the CMSSM model. The high scale 
parameters of the model were chosen to be m = 1460 GeV, = 600 GeV, 

tan f3 = 10, A Q = and \i > 0. This point lies just outside the exclusion limits 



set by the ATLAS experiment in 46 . The corresponding weak scale parameters 
and decay modes were calculated using ISA JET 7.80 (47). The Herwig++ im- 
plementation of the MSSM model was used to generate LHC pp collisions at a 
centre-of-mass energy of y/s — 8 TeV. In the following sections the impact of the 
POWHEG correction on a number of archetypal decays is presented. 



4.2.1 u L ^>-uXi 

Events were generated in which and its associated anti-particle were produced 
and then decayed via the mode ul — > v>Xi- The final-state particles from each 
decay, excluding the xl, were clustered into jets using the FastJet implementation 
of the anti-fc-r algorithm with the energy recombination scheme and R = 0.4. 
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Fig. 2: Comparison of distributions generated with and without the POWHEG 
correction for the decay ul — > u£i i n the CMSSM model with 
m = 1460 GeV, m l/2 = 600 GeV, tan (3 = 10, A = and fi > at 
the LHC with yfs = 8 TeV. The left-hand plot shows the logarithm of the 
jet measure y 2 \ and the right-hand plot the of the the second hardest 
jet in the decay, both in the rest frame of the u^. 

Jets with pt < 20 GeV or |^| > 2.8 were discarded. The jet measure y%i, as 



defined in 48 , was calculated in the rest frame of the decaying ul and the 
corresponding distribution of log y 2 \ is shown in the left-hand plot in Fig. [2j The 
right-hand plot in Fig. [2] shows the transverse momentum distribution of the 
second hardest jet in the decay. In both plots the blue dashed line corresponds to 
the distribution generated using the LO matrix element combined with the parton 
shower while the black solid line shows the result with the POWHEG correction 
to the decay applied. Including the POWHEG correction has the effect of shifting 
the distribution in the left-hand plot of Fig. [2] to lower values of the jet measure. 
This happens because the parton shower in Herwig+-|- tends to over-populate the 
hard region of phase space. With the POWHEG correction applied, the hardest 
emission in the shower is now correctly generated using the real-emission matrix 
element and hence the over population is reduced. The same effect can be seen 
in the px2 distribution in the right hand plot of Fig. [2j 

4.2.2 ul —> ug 

The scenario where ul decays through the mode Ul — >■ ug was also studied. 
To focus on the impact of the POWHEG correction to this decay the gluinos 
produced were not allowed to decay further. Events were generated at the parton 
level to accommodate the production of stable gluinos and were analysed as in 



Sect. 4.2.1 Distributions of log y^ 2 arid pT3 in the rest frame of the ul were 



plotted and are shown in the left and right hand plots of Fig. [3] respectively. 
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Fig. 3: Comparison of distributions generated with and without the POWHEG 
correction applied to the decay ul — >■ ug for stable gluinos. Results are 
for the CMSSM model with m = 1460 GeV, m 1/2 = 600 GeV, tan f3 = 10, 
Aq = and \i > for LHC events with yfs = 8TeV. The left-hand plot 
shows the logarithm of the jet measure y^ 2 and the right-hand plot the px 
of the the third hardest jet in the decay, both in the rest frame of the ul- 



The blue dashed line shows the purely LO result and the black solid line is the 



result including the POWHEG correction. As in Sect. |4.2.1[ we find that the 
POWHEG correction reduces the pt of the third hardest jet and shifts the jet 
measure distribution to lower values of log y^ 2 - 



4.2.3 b 2 -»■ W~ t x 

Finally, events were simulated where the heavier superpartner of the bottom 
quark, b 2 , was pair produced and decayed via b 2 — > W~ t\. Parton level events 
were generated with stable t\ and leptonically decaying W bosons. The final state 
particles were analysed as in Sect. |4.2.1 with the decay products of the W bosons 
omitted from the jet clustering. Distributions of logy 2 i an d Pt2 were plotted 
and are shown in Fig. [4} As in Sect. |4.2.1| and |4.2.2| we see that the POWHEG 
correction softens the log y 2 i and pr 2 distributions. 



5 Conclusions 

In this work we used the real-emission matrix element to generate hard QCD 
radiation in a range of particle decays in the Herwig++ event generator. This 
method is particularly relevant to new physics searches based on the decays of 
heavy new particles produced in LHC collisions. The POWHEG corrections to 
these decays can change the shapes of certain experimental observables, thus 
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Fig. 4: Comparison of distributions generated with and without the POWHEG 
correction applied to the decay b 2 — > W~ ti for stable t\ and 
leptonically decaying W bosons. Results are for the CMSSM model with 
m = 1460 GeV, m 1/2 = 600 GeV, tan (3 = 10, A = and \i > at the 
LHC with a/s = 8 TeV. The left-hand plot shows the logarithm of the jet 
measure y 2 i and the right-hand plot the pt of the the second hardest jet 
in the decay, both in the rest frame of the b 2 - 

altering the number of signal events passing selection cuts and modifying the ex- 
clusion bounds that can be set on the masses of the new particles. This correction 
will be available in a future version of Herwig++ . 

The technique was tested by plotting the invariant mass distribution of dijets 
produced in the decay of the lightest RS graviton, G — > gg and G — > qq. Applying 
the POWHEG correction was found to have a considerable impact on the height 
of the distribution in the dijet mass peak. The percentage of events passing 
selection cuts in the mass range 2.1 TeV < rrijj < 2.3 TeV dropped from 2.4% 
to 1.5% when the correction was applied. This corresponds to a decrease by a 
factor of approximately 1.6. This is a consequence of the dijet invariant mass not 
including the hardest emission in the shower that will carry a significant fraction 
of the graviton's momentum when it is simulated using the real-emission matrix 
element. This illustrates the importance of including higher order corrections 
when optimising experimental searches. 

The POWHEG correction was also tested for a range of decays in the CMSSM 
model. Here it was found that the correction reduced the transverse momentum 
of the hardest emission in the decay. This leads to a softening of the pt and jet 
measure distributions. 

In this work we have used the POWHEG formalism to improve the simulation 
of hard radiation in particle decays and studied the resulting effect on experimen- 
tal observables. However, hard radiation in the initial-state parton shower can 
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also have a significant impact on these distributions. Hence, in order to achieve 
accurate simulation of hard radiation in BSM processes we must also include ef- 
fects from the initial-state shower. Using the POWHEG formalism to improve the 
simulation of the hardest initial-state emission in the shower will be the subject 
of future work. 
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